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We report numerical and analytical studies of the reversal of the magnetic moment of a single- 
domain magnetic particle by a circularly polarized ac field of time-dependent frequency For the 
time-linear frequency sweep, the phase diagrams are computed that illustrate the dependence of the 
reversal on the frequency sweep rate v, the amplitude of the ac field h, the magnetic anisotropy field 
d, and the damping parameter a. It is shown that the most efficient magnetization reversal requires 
a non-linear time dependence of the frequency, cj(t), for which an exact analytical formula is derived 
with account of damping. The necessary condition of the reversal is h > ad. Implementation of a 
small-scale magnetization reversal is proposed in which a nanomagnet is electromagnetically coupled 
to two weak superconducting links controlled by the voltage. Dynamics of such a system is analyzed 
with account of the back effect of the magnet on the superconducting links. 

PACS numbers: 75.60. Jk, 84.40.-x, 75.50.Tt, 85.25.Cp 



I. INTRODUCTION 

In recent years a significant effort has been made 
to achieve magnetization reversal in nanostructures, as- 
sisted by the low amplitude ac field in the radio frequency 
range. The idea is rather simple. The dc magnetic field 
required to reverse the magnetization of a single-domain 
magnetic particle, the so-called anisotropy field, is typ- 
ically in the range 0.01-0. IT. The field of this strength 
at the location of the particle is not easy to develop fast. 
The ac magnetic field that one can typically develop in 
the radio-frequency range would be two orders of mag- 
nitude weaker. Applied in a resonant fashion, it could 
increase the amplitude of the precession of the magnetic 
moment, sometimes leading to its full reversal, in the 
same way as weak pushes of a pendulum at the frequency 
of its mechanical oscillation can flip the pendulum over 
the top. However, the study of both problem shows a lack 
of robust reversal. In many cases the attempted reversal 
results in a chaotic behavior that may be interesting on 
its own. 

Magnetization reversal by ultrashot magnetic field 
pulses produced by a high-energy electron beam has been 
studied by Back et al. 1 in perpendicularly magnetized 
CoPt films. Schumacher et al 2 studied phase coherent 
precessional magnetization reversal in spin valves by a 
pulse of the transverse field of a few hundred picosecond 
duration produced by the electric current. 

Later, significant number of experiments focused on 
microwave-assisted reversal in smaller structures and 
individual single-domain magnets with a strong static 
field applied to reduce the barrier. Thirion et aP at- 
tempted magnetization reversal in static fields below 
the anisotropy field, assisted by a linearly polarized mi- 
crowave field, in 20-nm-diameter Co particles placed on 
the bridge of a micro- SQUID. They were able to repro- 
duce the Stoner-Wohlfarth astroicP and study the de- 
pendence of the reversal on the frequency and duration 



of the ac pulse. Enhancement of the magnetization re- 
versal by microwave magnetic fields in nanometer Co 
strips has been demonstrated by Grollier et al. 6 Nem- 
bach et aP and Nozaki et aP used magnetic force mi- 
croscopy to measure microwave assisted magnetization 
reversal in individual submicron Co and permalloy par- 
ticles. Microwave-assisted magnetization switching in 
permalloy tunnel junctions has been demonstrated by 
Moriyama et al 9 . Podbielski et al studied magnetiza- 
tion reversal in microscopic permalloy rings at GHz fre- 
quency. They observed non-linear spin dynamics and 
obtained experimental phase diagram of the reversal 
as function of microwave frequency and power Using 
time-resolved magneto-optic Kerr microscopy, Wolters- 
dorf and BaclP^ detected enhancement of magnetization 
switching in single-domain permalloy elements subjected 
to the resonant microwave field. Microwave-assisted mag- 
netization reversal in single-domain permalloy nanoele- 
ments has been studied by Nembach et alP^ Wang et 
aP^ have investigated experimentally the competition 
between damping and pumping for microwave-assisted 
magnetization reversal in FeCo thin films. 

Theoretical research in this area mostly focused on the 
magnetization reversal assisted by the ac field of con- 
stant frequency!^ Non- linear magnetization dynamics in- 
duced by such a field that results in a c haotic behavior 
has been studied by Bertotti et al! 16 * 17 ! Denisov et aP^ 
addressed magnetization of nanoparticles in a rotating 
magnetic field. Synchronization and chaos induced in 
the damped dynamics of a single-domain particle by the 
ac field of constant frequency has been investigated by 
Sun and Wang. 5 Nonlinear-dynamical-system approach 
to the microwave-assisted magnetization dynamics was 
reviewed by Bertotti et alP^ Micromagnetic modelling 
of microwave-assisted magnetic recording was performed 
by Wang et al. 12 Constant frequency microwave switch- 
ing magnetic grains coupled by exchange interaction has 
been investigated by Igarashi et alP^ Okamoto et al ad- 
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dressed stability of the magnetization switching by lin- 
early and circularly polarized waves. Magnetization re- 
versal in a resonant cavity has been studied by Yukalov 
and Yukalova!^ 

Fewer number of theoretical papers have considered 
dynamics of the magnetization of a nanoparticle gener- 
ated by the ac magnetic field of variable frequency. May- 
ergoltz et al^ developed the inverse problem approach 
to the precessional switching of the magnetization by a 
linearly polarized pulse of the magnetic field. Rivkin and 
KettersorP^ obtained the optimal time dependence of the 
microwave frequency in the absence of damping, as well 
as the condition of the reversal in the presence of damp- 
ing. Magnetization reversal by a linearly polarized ac 
field of frequency that depends linearly on time has been 
studied by two of the authors.^ Barros et aP^" developed 
an optimization method in which the energy consump- 
tion needed for reversal is minimized with respect to the 
time dependence of the amplitude and frequency of mi- 
crowaves. 

A few general points need to be made before address- 
ing the problem of the reversal of the magnetization by 
the microwaves. Firstly, a robust magnetization reversal 
can be effectively achieved only with a circularly polar- 
ized ac field. Indeed, photons of circular polarization 
have a definite orientation of their spin projection, while 
photons with linear polarization are in a superposition 
of spin states. Consequently, photons with the right cir- 
cular polarization, when absorbed by the magnet, drive 
the magnetization in one direction towards the reversal, 
while linearly polarized photons can be both absorbed 
and emitted and, therefore, do not necessarily reverse 
the magnetization. Secondly, the photons are effectively 
absorbed only when they are in resonance with the spin 
levels. The latter are not equidistant on the magnetic 
quantum number, that is, on the projection of the mag- 
netic moment on the direction of the effective field. Thus, 
as the magnetic moment reverses, the photon frequency 
that can be resonantly absorbed by the magnet changes 
with time, so that the frequency of the microwave field 
has to be adjusted. Damping of the precession adds an- 
other dimension to this problem as the power of the ac 
field that is pumped into the magnet should exceed the 
rate of energy dissipation. Analysis shows that circularly 
polarized small-amplitude ac field of a time-dependent 
frequency that follows the condition of the resonance is 
sufficient for achieving magnetization reversal. The case 
of a zero static field is of the highest practical importance. 

The typical wavelength of microwaves that are in res- 
onance with the precession of the magnetic moment is 
in the centimeter range. Thus, one of the challenges for 
potential applications of the microwave- assisted magne- 
tization reversal for, e.g., computer technology, consists 
of the generation of a circularly polarized ac field of suf- 
ficient amplitude at the position of a nanoscale single- 
domain particle. In Ref. [25] a suggestion has been made 
to use the ac field generated by a superconducting weak 
link. If one is not turned off by the necessity to go 



to lower temperatures (which is probably, inevitable for 
magnetic memory of ultra-high density), the advantage 
of this method would be the possibility to control the 
time dependence of the frequency by voltage across the 
link. Interaction between a nanomagnet and a Josephson 
junction has been subject of intensive research. Micro- 
SQUID setup has been used by Jamet et al to ob serve 
switching of the magnetization in a 3nm Co clusteiP^S 
see also review 28 . Ferromagnetic resonance in permalloy 
films grown on Nb substrate has been studied by Bell et 
alP Petcovic et al 30 investigated experimentally the di- 
rect dynamical coupling of spin modes and a supercurrent 
in a ferromagnetic junction, following theoretical study 
of this system by HouzetP^ Current-phase relation in a 
Josephson junction coupled with a magnetic dot has been 
investigated theoretically by SamokhvalovP^ Most of the 
research in this area focused on the proximity effect^M^ 
rather than on electromagnetic interaction. 

In this paper we study magnetization dynamics of a 
single-domain uniaxial magnetic particle in zero static 
field, induced by a circularly polarized ac field of con- 
stant amplitude but variable frequency. The model is 
formulated in Section |TTJ General properties of the mag- 
netization reversal are studied in Section IIIII Numeri- 
cal results for the time-linear frequency sweep are pre- 
sented in Section [IV| where the phase diagrams are com- 
puted for the dependence of the magnetization switching 
on parameters. They are the frequency sweep rate, the 
amplitude of the ac field, the magnetic anisotropy field, 
and the damping parameter. Analytical results for the 
time-linear sweep, that are generally in good agreement 
with numerical results, are given in Section VA In Sec- 



tion VB we obtain the exact analytical solution for non- 
linear time dependence of the frequency that provides 
the fastest magnetization reversal. The model in which 
circularly polarized ac field is generated by two supercon- 
ducting weak links is studied in Section [VI] with account 
of the back effect of the dynamics of the magnetic mo- 
ment on the links. Our conclusions and suggestions for 



experiment are presented in Section VII 



II. THE MODEL 

The energy of a single-domain magnetic particle with 
an uniaxial anisotropy in a circularly-polarized ac field 
has the form 

H = -KVMl - VM x hcos<S>(t) - VM y hsm$(t). (1) 

Here K is the magnetic anisotropy constant, V is parti- 
cle's volume, M is the magnetization, h is the amplitude 
of the ac field, and is the phase related to the time- 
dependent frequency as 



$(t) = u(t). 



(2) 



One of the cases we consider is that of the frequency 
linearly changing with time, 



Oj(t) 



-vt, 



(3) 
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where = —vt 2 /2. The other case that will be studied 
here is a non-linear time dependence of the frequency that 
provides the fastest magnetization reversal. 

It is convenient to recast the problem in terms of a 
classical spin s = M/M s , |s| = 1, where M s is the sat- 
uration magnetization. The Landau-Lifshitz equation of 
motion for this spin has the form 



s = 7 [s x H eff ] - [s x [s x H eff ] 



(4) 



where 7 is the gyromagnetic ratio, a is dimensionless 
damping coefficient and 



H 



eff 



1 on 

VdM 



2ds z e z + he x cos + he y sin &(t) 

(5) 

is the effective field. Here d = H a = KM S is the 
anisotropy field. In the initial state the spin points in 
the negative- z direction, s(— 00) = — e z . 

Further it is convenient to switch to the coordinate 
frame rotating around the z axis together with the mag- 
netic field, so that in this frame the magnetic field is 
static. As the result, in the new frame the spin acquires a 
rotation opposite to that of the ac field in the initial (lab- 
oratory) frame. Thus in the rotating frame the Landau- 
Lifshitz has the form 

s = [s x ( 7 H eff + «(*))] - a 7 [s x [s x H eff ]] , (6) 

where 

H e ff = 2ds z e z + he x , ft(t) = u(t)e z . (7) 



III. GENERAL PROPERTIES OF THE 
MAGNETIZATION REVERSAL 

With the sign choice in Eq. (|3), the ac field at negative 
times is precessing in the same direction as the magnetic 
moment, thus it excites magnetic resonance and may 
cause magnetization reversal. In the ideal case, as we 
will see below, the resonance condition holds during the 
whole reversal. After the magnetic moment overcomes 
the barrier, s z > 0, it changes its precession direction, 
and so does the ac field. 

In the rotating frame, the field uj(t)/j sweeping at a 
linear rate makes the problem resembling that of the 
Landau-Zener (LZ) effect that can be formulated in terms 
of the evolution of a classical spin described by a Larmor 
equation. There are three modifications, however: (i) 
uniaxial anisotropy and (ii) damping added to the model 
and (iii) sweeping uj(t) in the negative direction. Because 
of the latter, the initial state of the spin in the rotating 
frame is the high-energy state with s z opposite to uj, see 
Eq. ([3]). To the contrast, in the regular LZ effect the 
initial spin state is the low-energy state. In the absence 
of anisotropy and damping, the initial orientation of the 
spin and the sweep direction do not matter. However, 
in the general case the situation does depend on these 
factors. 



In particular, in the absence of damping one can mul- 
tiply 7h e ff + by —1 that only makes the spin precess 
in the opposite direction but does not affect its reversal. 
The resulting model is a model with a positive sweep, 
such as the regular LZ problem, while the anisotropy be- 
comes easy-plane, d < 0. It was showri^ that in this 
case for a sweep slow enough the system adiabatically 
follows the time-dependent lowest-energy state and the 
spin switching is very efficient. In our original model 
(with a = 0) the magnetization reversal is similar. Only 
instead of adiabatically following the lowest-energy state, 
the spin adiabatically follows the highest-energy state, in 
which it was at the beginning. 

This adiabatic solution corresponds to the maximum 
of the energy in the rotating frame 



H/{VM S ) = -dsl - s x h - (w/7) s z 



(8) 



The maximal energy corresponds to s y = 0. Using s x = 
— \J\ — s 2 z (opposite to the transverse field) and requiring 
dl~L/ds z = 0, one obtains the equation 



2ds z + hs z /y/l-s 2 z + uj/j = 



(9) 



for the energy maximum. Since in practical conditions 
h <C d, an approximate solution of this 4th-order alge- 
braic equation for the adiabatic spin value reads 



-1, uj > 2>yd 

s z = <( -u/(2jd), \uj\ < 2jd 



(10) 



1, 



uj < -2jd. 



Note that this solution is independent of h. Nonzero 
values of h cause rounding at the borders of the cen- 
tral region |cj| < 2jd where the reversal occurs. In the 
laboratory frame, the spin is precessing during adiabatic 
reversal being phase-locked to the ac field. 

For a = the magnetization reversal can be achieved 
for whatever small ac field h. In the case of a nonzero 
damping, there is a dissipative torque acting towards the 
energy minima, and the ac field h has to exceed a thresh- 
old value to overcome this torque. Below we will see that 
the magnetization reversal requires 



h > ad 



(ii) 



that is much easier to fulfill than h > 2d in the case of 
a static field. As the torque due to the transverse field 
is maximal when the magnetic moment is perpendicular 
to it, in the presence of damping the magnetic moment 
goes out of the x-z plane during the reversal. 



IV. NUMERICAL - MAGNETIZATION 
REVERSAL BY THE TIME-LINEAR 
FREQUENCY SWEEP 

A. Time dependences of reversing magnetization 

The results of numerical solution of Eq. (|6| in the un- 
damped case a = for a small frequency sweep rate are 
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Figure 2: Non-adiabatic magnetization reversal at zero damp- 
ing. 



shown in Fig. [T] Magnetization reversal in this case is 
almost adiabatic and s z (t) is well described by Eq. (10) 



with rounding at the borders of the reversal interval due 
to a small value of h/d. The reversal is practically con- 
fined to the z-x plane and s y is small. Numerical results 
for a faster sweep rate are shown in Fig. [2] Here there 
is still magnetization reversal but it is not adiabatic and 
the final value of s z is smaller than one. Because of this, 
the magnetic moment is precessing around the z axis, as 
manifested by s x and s y . During reversal the magneti- 
zation is substantially deviating from the z-x plane. For 
larger sweep rates the reversal quickly becomes impossi- 
ble. 

Fig. [3] shows that in the damped case the magnetic 
moment substantially deviates from the z-x plane. Still, 
overall the reversal in this case is close to adiabatic. In- 
creasing the sweep rate leads to a non-adiabatic regime 
shown in Fig. [4j Here transverse spin components are os- 
cillating and the dependence of s z is jagged. This shows 
that, in the laboratory frame, the phase locking between 
the magnetic moment and the ac field is about to break. 
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Figure 4: Non-adiabatic magnetization reversal for a — 0.02. 



In spite of all this, there is a complete reversal because 
the damping finally brings the magnetic moment to the 
bottom of the potential well (c.f. Fig. |2|. For a faster 
sweep the reversal disappears and the magnetic moment 
lands in the initial well, s z = — 1. In the case of a slow 
sweep shown in Fig. [5] an instability can develop that 
leads to the breakdown of the phase locking and to faster 
relaxation of the magnetic moment towards one of the 
two potential wells. The final value of s z (1 or -1) be- 
haves irregularly vs sweep rate. This regime is not in- 
teresting for applications aimed at achieving as fast as 
possible reversal. 



B. 



Phase diagram of the magnetization reversal by 
the time-linear frequency sweep 



Dependence of the final value of s z on the amplitude 
of the ac field h and frequency sweep rate v defines the 
"phase diagram" of the magnetization reversal. In the un- 
damped case the numerically calculated phase diagram is 
shown in Fig. [6j The final s z is color-coded: black corre- 
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Figure 7: Phase diagram of t he magnetization reversal for 
a = 0.02. Yellow line: Eq. p6). 
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Figure 6: Phase diagram of the magnetization reversal in the 
undamped case. 



sponds to s z = — 1 (non-reversal) and red corresponds to 
s z = 1 (reversal). One can see that the reversal requires 
h sufficiently large and v sufficiently small. The curva- 
ture of the phase boundary at small h and v suggests a 
fractional power. Careful examenation of this region of 
the phase diagram shows that the reversal condition has 
the form 



2j 2 d 2 



< c 



-) 

d 



4/3 



1.6 



(12) 



at h/d <C 1 and a = 0. 

Phase diagram of the magnetization reversal in the 
damped case a = 0.02 is shown in Fig. [7] It is similar to 
Fig. [6] but there is a threshold for the magnetization re- 
versal on h and the phase-boundary line goes linearly at 
small v. Computations for different values of a suggest 
that the reversal requires h/d > a. 

One can compute other types of phase diagrams for 
the magnetization reversal that show a compact rever- 
sal region and the whole boundary line. The most 
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Figure 8: Efficiency-type phase diagram of the magnetization 
reversal for a — 0.1. 



useful of these phase diagrams uses the parameters 
(ad/h, av/^h 2 )). Indeed, the area of the magnetiza- 
tion reversal is the compact region < ad/h < 1 and 
v/h 2 is inversely proportional to the energy of the ac 
field injected during the time of the reversal by the lin- 
ear frequency sweep 



^(linear) _ ^7^ 



(13) 



The maximum of v/h 2 corresponds to the minimal in- 
jected energy and thus to the maximal efficiency of the 
reversal. Figs. [5] and [9] show that the maximal effi- 
ciency of the time-linear frequency sweep corresponds to 
ad/h « 0.5. Also in these figures one can see that there 
is no reversal if the sweep rate is too low, especially for 
low ac fields on the right side. 
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Figure 9: Efficiency- type phase diagram of the magnetization 
reversal for a = 0.01. 



V. ANALYTICAL 

Analytical investigation of the magnetization reversal 
is more convenient in spherical coordinates 



s z = cos 



sin U cos (/?, s y = sin U sin 



(14) 

After neglecting the ac field in the dissipation term, Eq. 
(|6| becomes 

= sin cp — ajd sin 20 (15) 
<£> = — 27<icos — uj (t) + 7/1 cos 99 cot 0. (16) 

A. Linear frequency sweep 

In the case of a linear frequency sweep, Eq. ([3|, one 
can rewrite the equation of motion for the spin in terms 
of the dimensionless time variable 

r = vt/(2jd). (17) 

The resulting equation of motion has the form 

dO/dr = b sin cp — aa sin 20 (18) 
dcp/dr = — 2a (cosO — r) + b cos (p cot 0, (19) 



where 



2 7 2 d 2 , 2 7 2 d/i 
a = , 6 = 



(20) 



characterise the sweep rate. Another important parame- 
ter is 



A — ad/h. 



(21) 



Since a is a large parameter, phase locking of the mag- 
netic moment to the ac field and thus efficient reversal 



requires cos# = r in the reversal region |r| < 1. If cos# 
only slightly deviates from this form, this causes strong 
oscillations of cp and thus the breakdown of the phase 
locking. Setting cos# = r, from Eq. (18) one obtains the 



phase-locking condition for cp in the form 



sin cp 



1 dO 

bah 



Asm 20. 



(22) 



The term on the left of this formula is proportional to 
the torque acting on the spin from the ac field. This 
torque has to ensure temporal change of (i.e., reversal) 
and compensate for the dissipative torque that is acting 
toward potential wells. One can see that damping ham- 
pers climbing the barrier by the magnetic moment. The 
maximal damping torque is realized at = 37r/4, where 
sin 2$ = —1. Since the reversal implies dO/dr < 0, it is 
clear that for A > 1 the ac torque cannot overcome the 
damping torque. Thus, the necessary condition for the 
magnetization reversal is 



A<1, 



(23) 



while the more restricting sufficient condition requires 
that the right-hand side of Eq. ( 22 ) does not drop below 



-1 for all r. The latter requires the frequency sweep rate 
to be not too fast. Using cos# = r, one can rewrite this 
condition in the form 



max/(r) < 1, f(r) = —2Ar\/ 1 - r 2 



bVT 



10 



(24) 
di- 



Because of the inertial term, f(r) shown in Fig. 
verges at the borders of the reversal interval. This is, 
however, an artefact of neglecting the rounding of the de- 
pendence s z (t) at r = ±1 because of the finite value of 
h. When this effect is taken into account, there are max- 
ima around r = ±1 instead of divergences. Thus, spin 
reversal can break down either because of the inertial ef- 
fect near r = — 1 or because of the effect of dissipation 
near = 37r/4, i.e., r = —l/y/2, depending on which one 
occurs at a smaller sweep rate. 
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Let us first consider the dissipative breakdown of the 
magnetization reversal for A slightly below 1 that hap- 
pens at a small sweep rate. In this case 1/b oc v is small 
and the second term in f(r) in Eq. ( [24] ) is a perturba- 
tion. Thus the value of this term can be taken at the 
unperturbed dissipative maximum r = —l/y/2. Using 
— 2rVl — t 2 = 1 and y/l — r 2 = 1/V% one obtains the 
reversal condition 



A + 4 <1. 
o 

This can be rewritten in real units as 



< 



27 2 <f 2 ^2 \ d 



(25) 



(26) 



and it is in a reasonable agreement with the numerical 
results in Fig. [7| Although this expression formally sur- 
vives in the dissipationless limit a — >• 0, it becomes inap- 
plicable in this limit. Here the breakdown of spin reversal 
is due to the inertial effect. 

To investigate the latter, one needs a more accurate 
approximation for f(r) in Eq. (24) near r = — 1 that 



tranforms divergence into a maximum. This can be done 
by solving Eq. (|9) although it is difficult to do it ana- 
lytically in general. Instead, since the maximum should 
be close to r = — 1, we can solve this equation exactly at 
r = — 1, which is much easier. A perturbative solution 
for h/d <C 1 yields 



ds/dr ^ 2/3, sin ^ (h/d) 



1/3 



and then one obtains 



dO 

dr sin dr 



1 dcosQ 2 ( d 



3 V ft 



1/3 



(27) 



(28) 



Replacing in Eq. ( 24 ) 1 / yl — r 2 by this result and using 
Eq. (20), one obtains the dissipationless reversal condi- 
tion 



v 3 /ft\ 
2 7 2 d 2 < 2 \d) 



4/3 



(29) 



in a reasonable agreement with the numerical result, Eq. 
The combined reversal condition obtained from Eqs. 



(26) and (12) is thus 



v . [ i (h \ 3 m 4/3 ~ 
2W <mm K7lU" a J'2Uj 



(30) 



Let us shortly discuss the stability of our spin-reversal 
solutions that, in the laboratory frame, is the stability of 
phase locking between the spin and the ac field at slow 
frequency sweep. Linearizing Eqs. (15) and (16) around 



the static solution (0, ip) at a fixed time, one obtains the 
deviation (50, Sep) oc e xt . For the orientations closer to 
the wells, 37r/4 < < 7r and < < 7r/4, one has A < 



and phase locking is stable. However, for the orientations 
closer to the barrier, 7r/4 < < 37r/4, one has A > and 
phase locking is unstable. Thus the barrier has to be 
crossed fast enough during reversal before the instability 
develops. Considering the process quasi-statically, one 
can write 



((5(9, Sep) ~ exp 



dt'\(t') 



Uo 



(31) 



and use the stability criterion J tQ dt\(t) < 1, where to is 
the time of entering the instability region and the top of 
the barrier is reached at t = 0. After some algebra one 
arrives at the stability criterium 



a 



< 



3V2 2 7 2 ^ 2 ' 



(32) 



A boundary of this kind is seen in Figs. [8] and [9] close to 
the bottom. 



B. Optimal frequency sweep 

The magnetization reversal can be optimized by apply- 
ing a time- nonlinear frequency sweep. Among all possible 
cases the so-called "optimal sweep" stands out as a rota- 
tion of the magnetic moment in one plane (in the rotating 
frame) with cp = — 7r/2, that is with the moment being 
always perpendicular to the ac field. It is easy to see that 
this provides the maximal torque on the magnetization 
during the reversal. With 



■(£) = -2-fdcosO 



(33) 



Eq. (16) self-consistently yields cp = 0. Then Eq. (15) 



takes the form 



= —7ft — a / ydsm29. 



(34) 



Integrating this equation with the initial condition 0(0) = 
7r, one obtains 



tan 



• sin (t) 



(35) 



where t = y/l — A 2 jht and A is defined by Eq. (21). 



After some algebra the expression for the optimal sweep 
can be transformed to the most convenient form: 



cos (t 



arcsm 



A) 



1 — A cos (2t + arccos A) 



(36) 



illustrated in Fig. 11 Together with Eq. (|33j) it gives 



a non-linear time dependence of the frequency of the ac 
field that provides the fastest reversal of the magnetic 
moment. This exact result, that generalizes the result of 
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£ 0.0 




Figure 11: Optimal magnetization reversal for different values 
of A = ad/h. 



Ref. f24l obtained in the absence of damping, must have 
important practical applications. 

In the dissipationless case, A — >> 0, the optimal mag- 
netization reversal is described by a pure cosine function 
that is a Rabi precession of the magnetic moment around 
the ac field. In the general case, the reversal is mainly 
due to the cos term in the numerator, whereas the de- 
nominator only affects the shape of the switching curve, 
making it non-sinusoidal in the presence of damping. 
Eq. 



(|36j) and Fig. 
during the time 
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describe the optimal reversal 



ad 

T 



(37) 



It is instructive to compare this time with the time of 
the magnetization reversal for the linear sweep, defined 
by Eq. (13) (notice that the total time of the process 
may be longer). In the dissipationless case, the maximal 
sweep speed is given by Eq. (12). For that speed one 



obtains the minimal reversal time 



^(linear) 



3jh \h 



1/3 



(38) 



0. 



that is longer than the time given by Eq. ( 37 ) with A 
In the dissipative case near A = 1, the maximal sweep 
rate for the linear sweep follows from Eq. (26). This 
yields 



t 



(linear) 



2^2 1 



(39) 



For |1 — A\ <C 1 this is also is much longer than the time 
given by Eq. (37). In the relevant region A ~ 1 Eqs. (37) 
and (39) are comparable. However, one has to keep in 



mind that the linear frequency sweep has to begin with a 
frequency beyond the resonance range, so that the actual 
reversal time of the linear sweep is somewhat longer than 
above. 

The total energy input of the ac power needed for the 
reversal satisfies 



E oc h 2 U 



(40) 



For the optimal sweep one has 

h h 2 

E oc 



Vl-A 2 



^h 2 - (adf 



(41) 



The minimum of this function, 2ad, is achieved at 

h = V2ad (42) 

(that is at A = ad/h = l/y/2). 
For the linear sweep one has 



^(linear) ^ /^(linear) = 



h 



h 2 



1 — A h — ad 
The minimum of this function, Aad, is achieved at 

h = 2ad 



(43) 



(44) 



(that is, at A = ad/h = 0.5). 

We see that the magnetization reversal by the opti- 
mal sweep requires both a smaller amplitude of the ac 
field and a smaller total energy input, as compared to 
the linear sweep. In both cases the injected energy at 
the maximal efficiency is proportional to a and thus the 
efficiency itself is inversely proportional to a. (The lat- 
ter was multiplied by a in Figs. [8] and [9] to make them 
approximately scale with a.) 



VI. REVERSAL OF THE MAGNETIZATION 
BY JOSEPHSON CURRENTS 

Pointed switching of the magnetization of a nanomag- 
net by the ac field of varying frequency may be achieved 
by coupling the magnet to a weak superconducting linkP^. 
Advantage of this method consists of the possibility to 
control the time dependence of the frequency by the volt- 
age across the link, V(t). As has been discussed in the In- 
troduction the most effective switching occurs when the 
ac field is circularly polarized. This requires two weak 
links shown in Fig. [l2| In addition to the previous prob- 
lem one should now take into account the back effect of 
the magnet on the weak links. As we shall see below, our 
results for the optimal sweep permit generalization that 
provides exact time dependences of voltages on the two 
links needed to obtain full reversal of the magnetization. 

Each superconducting weak link interacting with the 
magnet contributes the term 



£j = -Ej cos 



2tt 



dr • A(r, t) 



(45) 



to the total energy. Here Ej = hl c /(2e) is the Joseph- 
son energy of the link, with I c being the critical current. 
The argument of cosine is the gauge invariant phase that 
consists of two contributions. The first contribution sat- 
isfying S = 2eV(t)/h comes from the voltage across the 
link, while the second contribution is due to the vector 



1 vt/(2yd) 2 



Figure 12: Geometry used in the model: Nanomagnet makes 
the right angle with two parallel superconducting weak links 
of length 2L at a distance a from the magnet. 



Figure 13: Detrimental influence of the magnetization feed- 
back on the Josephson junction in the case of linear frequency 
sweep and zero damping. 



potential of the magnet integrated between the terminals 
of the link; <l>o = 2irHc/(2e) being the flux quantum. Fol- 
lowing the footsteps of Ref. [25j it is easy to show that 
interaction of the magnet with two Josephson junctions 
leads to a modified expression ([5| for the effective field, 

H eff = 2ds z e z + hj[sm(5 y - ks x )e x + sin(S x - ks y )e y ] 

(46) 

where hj = kEj/ (M S V) is the amplitude of the ac mag- 
netic field created by the junction at the position of the 
nanomagnet, 



2eV a 



x,y 



dt 



V Xi y are the voltages across the junctions, and 



(47) 



" a = 0.02 
- h Id = 0.03 
v/(2y 2 /z</) = 0.1 


/k =0.5^T 
fk = 




\k =0.1 

i i 



1 vt/(2yd) 2 



Figure 14: Influence of the magnetization feedback in the 
damped case with linear sweep. 



4ttM,V 



L 



a$ VL 2 + a 2 



(48) 



is a dimensionless spin- feedback coupling coefficient. 

For the time-linear voltage, dependence of the effective 
field in Eq. (46) on oscillating transverse spin compo- 
nents is detrimental for reversal. This happens because 
the oscillating additions to the otherwise time-smooth 
phases disturb phase locking between the magnetization 
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and the ac field and cause non-adiabaticity. Fig. 
shows that non-adiabaticity in the zero-damping case be- 
comes pronounced already for small values of the feed- 
back coefficient k. In the realistic damped case the nega- 
tive influence of finite k is even stronger. The instability 
of the phase locking for 7r/4 < < 37r/4 discussed at the 
end of Sec. V A| exponentially increases the mismatch 



between the directions of the ac field and the magneti- 
zation arizing because of the feedback. As a result, the 
magnetization randomly lands in one of the two wells, as 
shown in Fig. [14] 



In order to reduce the effect of the Josephson junctions 
on the magnet to the effect of a circularly polarized field, 
one can require that 

6 x -ks y = S(t), S y - ks x = S(t) + |, (49) 

where S(t) is a phase. This allows one to use the results 
of the previous section for the optimal sweep requiring 
S(t) = uj (t) = —2'ydcosQ. For such a sweep, s XjV (t) 
in the laboratory frame are precessing diring reversal as 
follows 



s x = sm0(t) cos S(t) 
s y = sm0(t) sm5(t) - 



sin(t) cos S(t) 



\f \ — A cos(2t + arccos A) 
sin(t) sin S(t) 



y/l — A cos(2t + arccos A) 



, (50) 



where A = ad/hj, t = \/l — A 2/ yhjt, and sm0(t) 



was 
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obtained by combining Eqs. (35) and (36). The time- 
dependent phase is given by 

S(t) = f Lj(t')dt' = 2d/kj f cos 0(P)dP . (51) 

Jo v 1 - A 2 Jo 



with uj(t) defined by Eqs. (33) and 8(t) given by (36). In 
accordance with Eq. (49), the optimal time dependence 



of the voltages across the junctions become 



V x (t) 



h 

Ye 



u(t)+ks y ] , V v (t) = —[uj(t)+ks x ]. (52) 

This dependence is shown in Fig. 15 The time depen- 
dence of s z in the figure is the same as for the optimal 
freqency sweep in Fig. [II] Oscillations of the voltages are 
due to the terms in Eq. (52) that depend on s XiV . They 



are weak as long as k is small. Oscillations disappear in 
the limit of k 0, making V XjV in that limit to follow the 
smooth time dependence of the optimal frequency sweep 
obtained in the previous section. 

Having practical applications in mind, it is interesting 
to test the stability of the reversal described by the above 
equations against high-frequency voltage noise. This can 
be done by writing 



VL,y{t) = V x ,y{t) 



e^F x 
e 



.»(*) 



(53) 



with F XjV being uniform random functions of time be- 
tween — 1 and +1 and e representing the relative strength 
of the noise. Quite remarkably, as is illustrated in Fig. 
[l6j the full magnetization reversal may occur even in the 
presence of a strong noise. At e = 1 this happens with 
more than 0.99 probability. With less than 0.01 prob- 
ability the magnetic moment bounces back to = —1 
before it reaches s z = 1. This can be traced to the fact 
that the high-frequency noise in most cases averages out 
in the phase 5 because the latter is proportional to the 
time integral of the voltage. 



VII. DISCUSSION AND CONCLUSIONS 

We have studied numerically and analytically a 
microwave- assisted reversal of the magnetic moment of a 
single-domain magnetic particle. We conclude from our 
studies that a circularly polarized ac field that has spe- 
cific time dependence of the frequency can be an effective 
tool for switching the magnetization. The corresponding 
physical mechanism consists of the resonant absorbtion 
of photons of the spin projection that ensures the consis- 
tent change in the projection of the magnetic moment. 
Emission of excitations by the magnetic moment inhibits 
this process. In the micromagnetic theory it is described 
by the phenomenological dimensionless small parameter 
a that can be independently measured. In single-domain 
particles this parameter is usually greater than in bulk 
materials and is typically of order 0.01 — 0.1. 37 The con- 
dition on the power of the ac field needed to overcome 
damping and reverse the magnetization is h > ad. For, 
e,g., the anisotropy field d of order 0.01T and a of order 
0.01, one obtains h of order 0.000 IT for the amplitude of 
the ac field, which is a reasonable value from the practical 
point of view. 

We have studied linear and nonlinear time dependence 
of the frequency of the ac field. It has been demonstrated 
that the linear case, uo = —vt, resembles the Landau- 
Zener problem. Magnetization reversal has been demon- 
strated numerically and the phase diagrams have been 
obtained that show the range of v, h, d, and a that 
provide the reversal. They show that for the reversal 
to occur, the frequency sweep must be sufficiently slow, 
but not too slow when the damping is finite. The linear 
case has also been studied analytically. Condition (30) 



has been obtained for the upper bound on the frequency 
sweep rate. For the values of the parameters used above, 
that upper bound is in the ballpark of 10 T GHz/s. The 
minimal reversal time for the time-linear sweep is of order 

( 7 h)-\ 
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We have also studied a time-nonlinear frequency sweep. 
Exact analytical solution for uo(t) that provides the 
fastest reversal has been obtained with account of damp- 
ing. It is given by equations (33) and (36). This finding 
may have important practical application. We call this 
sweep the optimal sweep. It has been demonstrated that, 
besides ensuring the fasted magnetization switch, it also 
pumps less energy into the system as compared to the 
linear sweep. In both cases the injected energy is pro- 
portional to a. 

Circularly polarized ac field can be generated by cou- 
pling a single-domain particle electromagnetically to two 
weak superconducting links whose phases are displaced 
by 7r/2 with respect to each other. One advantage of such 
a system is that the time dependence of the frequency of 
the ac field generated by the links can be controlled by 
voltage. This problem has been studied by us with ac- 



count of the back effect of the magnetic moment on the 
links. Magnetization reversal has been demonstrated nu- 
merically and analytical expressions have been derived 
for the time dependence of the voltages across the links 
that provide the fasted magnetization reversal. One re- 
markable property of this system is weak dependence of 
the reversal dynamics on the voltage noise. 
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